Draft version February 1, 2008 

Preprint typeset using lAT^^X style cmulateapj v. 04/03/99 



HYPERCRITICAL ADVECTION DOMINATED ACCRETION FLOW 
G. E. Brown, C.-H. Lee 

Department of Physics & Astronomy, State University of New York, Stony Brook, New York 11794, USA 

H.A. Bethe 

Floyd R. Newman Laboratory of Nuclear Studies, Cornell University, Ithaca, New York 14853, USA 

Draft version February 1, 2008 



o 
o 
o 

< 
(N 

> 
(N 
m 

O 

o\ 
o 



ABSTRACT 

In this note we study the accretion disc that arises in hypercritical accretion of M ~ 10* MEdd onto 
a neutron star while it is in common envelope evolution with a massive companion. Such a study was 
carried out by Chevalier (1996), who had earlier suggested that the neutron star would go into a black 
hole in common envelope evolution. In his later study he found that the accretion could possibly be held 
up by angular momentum. 

In order to raise the temperature high enough that the disc might cool by neutrino emission, Chevalier 
found a small value of the a-parameter, where the kinematic coefficient of shear viscosity is v — aCsH, 
with Cs the velocity of sound and H the disc height; namely, a ~ 10^^ was necessary for gas pressure 
to dominate. He also considered results with higher values of a, pointing out that radiation pressure 
would then predominate. With these larger a's, the temperatures of the accreting material are much 
lower, < 0.35 MeV. The result is that neutrino cooling during the flow is negligible, satisfying very well 
the advection dominating conditions. 

The low temperature of the accreting material means that it cannot get rid of its energy rapidly by 
neutrino emission, so it piles up, pushing its way through the accretion disc. An accretion shock is 
formed, far beyond the neutron star, at a radius ^ 10* cm, much as in the earlier spherically symmetric 
calculation, but in rotation. Two-dimensional numerical simulation shows that an accretion disc is 
reformed inside of the accretion shock, allowing matter to accrete onto the neutron star with pressure 
high enough so that neutrinos can carry off the energy. 

Subject headings: accretion, accretion disks - binaries: close - stars: neutron - black hole physics - 
stars: evolution 



1. INTRODUCTION 

In one stage of the theory of evolution of compact ob- 
jects the neutron star from the explosion of the primary 
progenitor enters into common envelope evolution of the 
evolving secondary. Recently, developing ideas suggested 
by Chevalier (1993) and later by Brown (1995), Bethe 
& Brown (1998) introduced hypercritical accretion into 
a population synthesis of binary compact objects. They 
found that in the standard scenario of compact binary evo- 
lution the neutron star essentially always went into a low- 
mass black hole. 

In scrutiny of hypercritical accretion during common en- 
velope evolution. Chevalier (1996) discussed the effects of 
rotation. He studied accretion discs and, in particular, 
neutrino-cooled discs and advection dominated discs. He 
chiefly considered discs with viscosity parameter a 10~^. 
Such small a's were needed in order to reach high enough 
temperatures, in situations where gas pressure predomi- 
nates, so that the temperature was sufficient for neutrino 
emission. 

Calculations (Brandenburg et al. 1996; Torkelsson et 
al. 1996) suggest that the higher values of a ~ 0.1, which 
Shakura & Sunyaev (1973) suggested as arising from mag- 
netic turbulence are more appropriate. Chevalier (1996) 
notes that "Radiation pressure tends to dominate at higher 



ing that radiation pressure predominates to the extent that 
we can take the adiabatic index to be 7 = 4/3. We show 
the resulting hypercritical accretion to be essentially that 
of Houck & Chevalier (1991) and Brown (1995), although 
angular momentum was not discussed in these earlier pa- 
pers. 



2. ADVECTION 

We consider super-Eddington accretion, so electromag- 
netic radiation is negligible compared to the energy carried 
by the infalling matter. We shall show that, for reasonable 
parameters, the temperature remains sufficiently low that 
neutrino emission is also small. Therefore the accreting 
matter keeps nearly constant energy per unit mass. 

The accreting matter orbits around the neutron star. 
Because of collisions between the accreting particles, we 
may consider their orbits circular. From Narayan & Yi 
(1994), the orbital velocity is 



(1) 



values of a 
r = 10* 

In this note we discuss the situation for a 



for M = 1.4 Mq, M = 1 
cm, the crossover point is a ~ 10^ 



yr 



and 



0.05, find- 



with vk = \/ GM/r is the Keplerian velocity and M the 
mass of the neutron star. 

Because of viscosity, the velocity has a radial (inward) 
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component which is conventionally set equal to^ 
dr 3 



(2) 



We shall assmiie a — 0.05. 

The potential energy per unit mass is —GM/r, and from 
eq. P) the kinetic energy from the bulk motion of disc is 
^ i?r = GM/7r. Since we assume negligible heat loss, 
the heat energy in the local frame rotating with the disc 
is 



3 6GM 
e = - X , 

4 7 r ' 



(3) 



where the factor 6/7 comes from the sum of the local po- 
tential and kinetic energies in the disc, and 1/3 of e hav- 
ing gone into PV work (Enthalpy rather than energy is 
conserved). The density is related to the rate of energy 
accretion by 



M ^ 2nr 2Hp (-f) 



(4) 



where H is the height of the accretion disc. An explicit 
solution of the accretion problem, by Narayan & Yi (1994) 



ir which we adopt. Then, with cq. 



finds H = 



M = 2'Kr [^—avKj P- 

Using eq. (|l|), and assuming a constant 

p cx r~^/^. 
The thermal energy per unit volume is 



pe = J cx r / , 



(5) 



(6) 



(7) 



which also gives the r-dependence of the pressure, p — 
pe/3. Let us assume the mass of a neutron star to be 
1.5 Alp, and its radius 10 km,^ then 



= 1.2 X lO^^cm^i s"^ 



(8) 



Let us further assume M = 1 Mq yr ^ = 6.3 x 10'^''^ g s ^, 
and a = 0.05, 



pt = 4.5 X 10^'^erg cm 



(9) 



The energy density of electromagnetic radiation plus elec- 
trons and positrons is 

pe = ^aT* = ^(1.37) x lO^^ T^^^ ergs cm-^. (10) 



Equating the two expressions 



T = 0.33 MeV. 



(11) 



This is only a moderate temperature. The emission of 
neutrino pairs, according to Dicus (1972) and Brown & 
Weingartner (1995) is 



e„ = 1.0 X 10 

= 4.6 X 10^*^ ergs s"^ cm~^ 



(12) 



This is the energy emission at the surface of the neutron 
star. It can easily be seen that this quantity goes like 
j„-45/8^ Multiplying by the volume element Airr'^dr and 
integrating over r, we get for the total energy loss by neu- 
trinos. 

Wn^ J en4:nr^dr = 2.2 x 10^^ ergs s'^ (13) 

On the other hand, the energy advected by infalling matter 
is 

Me = 6.3 X 10^^ g s~i x 1.3 x 10^° erg g"^ 

= 8.2 X 10"^ erg s~^ (14) 



So the fraction of energy radiated in neutrinos is 



Me 



= 2.7 X 10" 



(15) 



This proves our contention at the beginning of this section 
that the energy lost during infall in neutrino emission is 
negligible. 

If the rate of infall is changed or a is altered, the fraction 
of neutrino radiation changes as 



(16) 



Only if a were reduced by a factor lO'^'^, to about 2 x 10"^, 
would neutrino radiation become appreciable. 

Whereas the above might appear to solve the accretion 
problem, it will not give a steady state solution because 
the accreted matter, which is at too low a temperature to 
cool sufficiently rapidly by neutrino radiation, cannot be 
fitted into the narrow gap between neutron star and in- 
ner disc. Thus, it pushes outwards through the accretion 
disc (Chevalier 1999). This was first seen in the numerical 
calculation of Zel'dovich et al. 1972. In the spherically 
symmetric case an accretion shock must be built into the 
flow. Freely falling gas outside of this shock goes into sub- 
sonic flow inside, with formation of an envelope at small 
radii hot enough so that neutrinos can be emitted to carry 
off the energy. Detailed calculation (Houck & Chevalier 
1991; Brown 1995) gives 



rsh ^ 2.6 X lO^cm 



M 



-0.37 



Mq yr 



(17) 



^Our a here is (2/3) of the ass of Shakura & Sunyaev (1973), who took the kinematic coefhceint of shear viscosity to be = ^OssCsH. Thus 
eq. is that of Narayan & Yi (1994), u = aCsH. Our treatment is equivalent to their advection dominated solution with their / = e' = e = 1. 

^For the neutron star with M = 1.5 Mq, the Schwarzschild radius is Rsch = 4.4 km, and the marginally stable orbit is Rms = 3Rsch, 
which is larger than the radius of neutron star. From Abramowicz et al. (1988), however, the inner disk can be extended to the marginally 
bound orbit Rmb = "^Rsch with hypercritical accretion. In this note, therefore, we assumed that the accretion disc extends to the surface of 
the neutron star. 
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Chevalier (1999) illustrates the necessity for the accretion 
shock simply. Starting from hydrostatic equilibrium in the 
gap region 



dP _ 
dr ^ 



GMT 



(18) 



with r initially set to R because the gap between neutron 
star and the innermost part of the accretion disc is narrow, 
and r = (1 — Rsch/T)~^ 1 one can integrate from point A, 
the inner edge of the HADAF, to point B, basically at the 
outer edge of the neutron star. Taking 



P = Kp^ 



(19) 



since the flow is adiabatic, and integrating from point A 
to point B near the neutron star surface, one finds that 



(-B - (-A 



7- 



1 



GMT 
~R2- 



{Ra - Rb) 



(20) 



Now from eq. (|) eA ^ GM/R. For e to increase by a fac- 
tor ^ 2 so that neutrinos could carry off the energy in 
going from point A to B would require {Ra — Rb) ^ R, 
violating the assumption that the gap is narrow. This ar- 
gument is general, not depending on the dimensionality of 
the description. 

The above discussion is at least roughly consistent with 
the results of Armitage & Livio (1999) who carry out 2D 
numerical simulation of the hypercritical accretion. They 
find that an accretion disc forms inside the accretion shock 
and with the usually considered values of the viscosity the 
disc is probably able to transfer mass inwards at the Bondi- 
Hoyle (hypercritical) rate. Pressures at the neutron star 
are high enough so that the accreted energy can be car- 
ried off by neutrinos. They warn, however, that the pres- 
sure of a disc also makes the formation of outflows or jets 
probable, which could reduce the accretion rates onto the 
neutron star. 

Our results up through eq. ( p^ ) are consistent with those 
of Minishige et al. (1997) for black hole disc accretion, 
once our accretion rate is scaled down to the ~ M/MEdd ^ 
10^ that they use. For a black hole accretion energy can 
be lost from the system because the energy can be ad- 
vected across the horizon. However, even in the case of a 
black hole, an accretion shock is developed in the numer- 
ical calculations of MacFadyen & Woosley (1999). Their 
rate of accretion is ~ 0.1 Mq s~^, ^ 10^ times ours. In 
their evolution of the accretion first the low angular mo- 
mentum material in the equator and material along the 
axes falls through the inner boundary onto the black hole. 
An accretion shock at 350 km with an interior cen- 
trifugally supported disc then forms. The MacFadyen & 
Woosley (1999) calculation deals with a different situation 
from ours; namely, advection of material that is originally 
He onto a black hole. Here nuclear dissociation energy 
as well as neutrino emission plays an important role in 
lowering the pressure. From the Minishige et al. (1997) 
calculations at an accretion rate of M lO^Af Edd into a 
black hole one can see the formation of the accretion disc, 
but cannot discern an accretion shock. 

^Near the black hole, there should be relativistic corrections to the 
however, we assumed that the ADAF solution can be extended to the 



From the calculation most relevant to our work, that of 
Armitage & Livio (1999) one sees that the case the flow has 
initial angular momentum, the main effect of the angular 
momentum in their 2D simulation is to set both accretion 
shock and accretion disc into rotation. The matter then 
being advected smoothly onto the neutron star. Thus, 
modulo the possibility of jets, the resulting accretion is 
very much like that in the one-dimensional calculations of 
Houck & ChevaUer (1991) and Brown (1995). 

It should be noted, however, that in the one-dimensional 
case the pressure drops off as p ex r~^, appropriate for 
a r = 4/3 accretion envelope (Chevalier 1989; Houck & 
Chevalier 1991; Brown & Weingartner 1994), whereas in 
the ADAF solution p cx r^^/^. The flatter profile in the 
latter case makes it more difficult to reach the energy den- 
sities necessary for neutrino emission. 

Perhaps the most important result is that the orbital 
motion of the accreting material is preserved. So the neu- 
tron star, in the natural course of events, acquires an an- 
gular momentum. 



J = AM X r^VN 



(21) 



where AM is the total mass accreted, r^v is the radius of 
the neutron star, and vn the corresponding orbital veloc- 
ity, 



r^VN 



-GMr 



7 



N 



1/2 



0.76 X 10^^ cm^ s~^ (22) 



an appreciable angular momentum. 

These consideration apply equally to accretion to a black 
hole. The thermal energy, in this case, is not radiated 
away, but et^c"^ is added to the mass of the black hole. 
If we assume that the black hole has a mass Mb — 2 Mq , 
its Schwarzschild radius is 



rs = 6km. 



(23) 



From Abramowicz et al. (1988), for the hypercritical ac- 
cretion, the inner disk can be extended to the marginally 
bound orbit Rmb — '^Rsch — 12 km."^ So the orbital 
velocity at the inner disc boundary is 



VB = 0.80 X 10^"cm s~\ 
and the specific angular momentum 

rBVB 0.96 X lO^^cm^ s'^ 



(24) 



(25) 



about the same as for a neutron star. Going over to a 
Kerr Geometry with decreased event horizon this will be 
lowered somewhat. Bethe and Brown (1998) have consid- 
ered a neutron star of 1.4 Mq in common envelope with 
a giant, and have found that the neutron star will accrete 
about AM = 1 Mq of the giant's mass and thereby be- 
come a black hole. This accreted material will have an 
angular momentum 



Jb = AMtbVb- 



(26) 



Assuming that the core of the giant does not merge with 
the black hole, this core will later become a supernova, 

ADAF (Advection Dominated Accretion Flow) solutions. In this note, 
marginally stable orbit. 



4 



Hypercritical Advection Dominated Accretion Flow 



yielding a second neutron star. Gravitational waves will 
be emitted, and consequently later the two compact ob- 
jects will merge. It is important that the BH has a big 
angular momentum Jb which may be extracted (Bland- 
ford & Znajek 1977; Lee, Wijers, & Brown 2000), and 
may possibly give a gamma ray pulse: As we have seen, 
no special mechanism is needed to give that angular mo- 
mentum to the BH, but it arises naturally as the black 
hole is formed. 

3. DISCUSSION 

Although there is considerable uncertainty in the a- 
parameter in the viscosity, we believe our value to be in a 
reasonable range. Torkelsson et al. (1996) have analyzed 
simulations of the magnetic turbulence different groups 
used to estimate the strength of the turbulent viscosity. 
Estimates of the a-parameter which range from 0.001 to 
0.7 in the Shakura & Sunyaev (1973) definition, were ver- 
ified using the same code for all simulations. The higher 
values of a are the result of an applied vertical (roughly 
equipartition) magnetic field. Without this field the typi- 
cal value is a ~ 0.005. 

Armitage & Livio (1999) have raised the possibility of 
jets removing the infalling material, so that it will not ac- 
crete onto the neutron star. Jets are prevalent in situations 
with accretion discs, but the mechanisms for driving them 
are poorly understood (Livio 1999). Suggestions have been 
made that accretion discs could generate a magnetically 
driven outflow, first by Blandford (1976) and by Lovelace 
(1976). 

If the magnetic energy above the disc is larger than the 
thermal and kinetic energy densities, the ionized outflow- 
ing material is forced to follow field lines, some of which 
form an angle with the perpendicular to the disc surface. 
The material may then be accelerated by the centrifugal 
force like a bead on a wire, as described by Livio (1999). 

It would seem reasonable for jets to be driven in the 
range of accretion rates MEdd to lO^'MEdd, the latter be- 
ing the lower limit at which neutrinos can carry off the 
energy sufficiently rapidly that hypercritical accretion can 
take place. In this accretion interval, only Msdd can be 
accreted, and since the surplus material must leave, it 
seems reasonable that magnetic forces can coUimate it in, 
at least, some cases. We believe that the jets in Cyg X-3 
may be an example. 

It is more difficult to envisage such jet propulsion with 
hypercritical accretion, especially at rates close to M ^ 
10^ A^Edd where the Alfven radius comes right down to 
the surface of the neutron star; i.e., where the ram pres- 
sure always exceeds the magnetic pressure, unless the lat- 
ter is increased by differential rotation, etc. Of course, as 
in Eggum et al. (1988) who calculated only slightly super 
Eddington accretion {M ~ 4 A^Edd) a jet may be formed 



up the center of the disc axis, where there is little infalling 
matter, but this would not drastically lower the accretion. 

In fact, in a similar case of hypercritical accretion we 
have some observational support for the main part of the 
matter being accreted. In the evolution of the transient 
sources (Brown, Lee, & Bethe 1999) the H envelope of an 
~ 25 Mq ZAMS is removed in common envelope evolution 
with an ^ 1 Mq main sequence star. Following this, the 
Fe core of the massive star goes into a black hole of mass 
~ 1.5 Mq and most of the He envelope must be accreted 
in order to obtain the final ~ 7 Mq black hole. In general 
the He envelope will be rotating. In fact, in the case of 
GRO J1655-40 arguments have been made (Zhang, Cui, & 
Chen 1997) that the black hole is presently rotating close 
to the Kerr limit although later studies (Sobczak et al. 
1999; Gruzinov 1999) lower the velocity to ~ 0.6 - 0.7c. 
Because the ~ Rq radius of the He envelope is ^ 10^ 
times larger than the neutron star radius, even with slow 
rotation of the envelope an accretion disc will be formed 
before the matter can be advected into the black hole. 
The central He density is p 10^ g cm~^, about the same 
as the H density in our work above in the inner disc, for 
M ~ lO^AfEdd- Thus, the difference is that the matter is 
He, rather than H, in the transient sources. Clearly jet for- 
mation does not strongly affect the amount of He accreted 
onto the black hole. 

It should also be realized that observed masses of neu- 
tron stars are not very far below the Brown & Bethe (1994) 
proposed upper mass limit of ~ 1.5 Mq. Thus, only a 
small fraction of the Bethe & Brown (1998) accretion of 
~ 1 Mq need be realized for the neutron star to go into a 
black hole. 

Since entropy gradients become negative, advection 
dominated flows are unstable to convection, as found by 
Eggum et al. (1988) and as worked out in considerable 
detail by Narayan & Yi (1995). Advection dominated ac- 
cretion is found to dominate in a thick belt about the 
equator, with somewhat reduced velocity, depending upon 
the polar angle. Narayan & Yi do find that jets can be 
driven in fairly wide angles about the poles, but in agree- 
ment with Eggum et al. (1988) find that very little mass 
is driven off in these jets. 

Eq. (|^) shows the energy density to go inversely with 
a, as a~^. Thus a much higher temperature is reached 
for small a's, ~ 10 times greater for a ~ 10~^ than for 
a ~ 0.1 (Chevalier 1996), the energy density going as T^. 
The rate of cooling goes as oc a"^/** and increases a fac- 
tor ^ 10^ for the small a's. In this regime of temperatures 
the disc is cooled very rapidly by neutrino emission. 

We would like to thank Roger A. Chevalier and Insu Yi 
for useful suggestions and advice. This work was partially 
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No. DE-FG02-88ER40388. 



APPENDIX 

APPENDIX: ADVECTION DOMINATED DISCS 
In this Appendix we give a more conventional derivation of our results. 

In the case that the disc is thin to neutrino emission, the steady state disc equation can be written as (Shakura & 
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Sunyaev 1993; Shapiro & Teukolsky 1983, section 14.5) 

E = 2Hp, (Al) 

M = 2'KrY.Vr = const (A2) 

e-Kr^Hap ^ MiGMrf^ (A3) 
2,MGM 

/ X 1/2 / 3 X 1/2 

where we use the notation of Chevaher (1996). Q Here S is the surface density of the disc. Our equation of state is 

P^^— (A6) 

where we neglect the gas pressure. We assume that the neutrino pair process (Dicus 1972) is the dominant source of 
neutrino coohng, using the Brown & Weingartner (1999) expression 

e„ = 1.0 X 10^5 (^^) ergs s"! cm-^. (A7) 



Eq. (A4) is however, not correct for the main part of our accretion disc with radiation dominated pressure, because the 
temperature of the material is not high enough for it to be cooled appreciably by neutrinos. In the absence of appreciable 
cooling, we can use the Narayan & Yi (1994) advection dominated accretion, with self similar solution: 



3a (GM\ 



1/2 



5 + 2e V T 

5 + 2eJ \ J 

(5 + 2ef/H'I 



1 , (A8) 



(AlO) 



2 ■ 



25/237rar3/2(GM)i/: 
21/2(5 + 2e)i/2M(GM)i/ 



2 



P~ ^ 1 ' 5/2 — (All) 



where we take e = 1 for the case of 7 = 4/3. 

Energy is produced by viscous heating at a rate 



where we use H = from eqs. (A5,A10 All). In order to use eq. (A7) we must obtain a temperature. For a — 0.05, 
M = 1.5 Mq, and M = 1 Mq yr'^ = 0.63 x 10^^ g s^^ we find 

ep^„^ = 0.28 X 10^" r^^ ergs cm"^ s^^. (A13) 

Using 



dt 7 

we can integrate this expression: 



dr 3 GM 



Eprod = / fpioddt = / — — dr 



3MGM _7/9 2 r 3MGM 



'^aVGM8Tr{H/r) 5 3^ [gm Sirr^H 



* In this note, we take the kinematic coefficient of shear viscosity v = acH, while u = %olcH is used in eq. (17) of Chevalier (1996). 
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where the last expression shows epiod to be (2/5) of the value that would be obtained from dimensional 
quantity {2/b)r /\vr\ can be interpreted as a dynamical time 



Td = 



2 r 

5 \Vr\ 



= 1.3 r. 



3/2 



sec 



and 



Cprod = Td Eprod = 3.6 X 10 " Tg " ' ergs cm 
We obtain a temperature of the advected material by relating T to the thermal energy density 



analysis. The 
(A16) 

(A17) 



eth 



3.8 X 10^^ T^^y ergs cm'^. 



From eq. (A17) we find I 



T^,y = 0.96 X 10 



-7 -5/2 



showing that T'* goes as 



so that 



-5/2 



The total cooling by eq. (A7) will be 
Se = Td in 

Se e„ 0.29 x 10" 



(A18) 



(A19) 



(A20) 



f^prod Cprod 



13/8 



Thus i5e/eprod ^ 1-6 x 10^^ for r — 10^ cm at the inner disc just outside the neutron star, in rough aj 
eq. (p^). In other words, temperatures of the advected material are so low that the neutrino cooling is ne§ 
the material will be deposited essentially adiabatically onto the neutron star, heating it. 



(A21) 

;reement with 
;ligible. Thus, 
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